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Abstract 

The contributions to heavy meson mass differences by the strong hyperfine interaction, the light 
quark masses and the electromagnetic interaction are obtained from the empirical values of the D, 
D*, B and B* masses by means of a mass formula based on the heavy quark mass expansion. The 
three different types of contributions are determined with significant accuracy to next to leading 
order in that expansion. 



1 



I. INTRODUCTION 



Hadron masses reveal key aspects of QCD. The light pseudoscalar meson masses reveal 
the spontaneous breaking of chiral symmetry as well as its explicit breaking by the light 
quark masses, and in a rather direct manner they permit to extract the ratios of these 
quark masses which are fundamental parameters of QCD. They also give access to the 
electromagnetic contributions to the masses, for instance, in the mass difference between 
charged and neutral pions. In the framework of chiral perturbation theory, a rather accurate 
understanding of the effects of light quark masses and electromagnetic corrections on the 
pseudoscalar octet has been achieved, and to a lesser extent in baryons as well. Heavy mesons 
represent another kind of hadronic system where one can arrive at a good determination 
of the various effects that determine their masses. Also, to a first degree of approximation 
in the heavy quark expansion, the strong hyperfine effects, which involve the heavy quark 
spin, serve to determine the ratio m c /mb, which is another fundamental input in QCD. 
In this work, we show that the current knowledge of the heavy meson masses allows for a 
quantitative determination of the different effects that contribute to the differences of heavy 
meson masses. The approach followed here is similar to the one given by Rosner and Wise 
and the improvement in the results is possible thanks to the better empirical accuracy 
in the heavy meson masses and also in the better knowledge of the heavy quark masses. 

In the limit of infinite heavy quark masses and of light quark SU(3) symmetry, heavy 
mesons fill multiplets of U(2Np) x {7(3) x SUj(2), where Np = 2 is the number of heavy 
quarks, SUj(2) is the rotation group associated with the light degrees of freedom of the 
heavy meson. In particular, the ground state mesons, namely D, D*, B and B*, fill the 
multiplet (4, 3, 2) of that symmetry. The symmetry is broken by several effects: 

• The finite masses of the heavy c and b quarks break U (4) x SUj(2) — > SUj(2) x U c (l) x 
{7,(1), where SUj(2) is the rotation group is associated with the spin of the meson. 

• The light quark masses break {7(3) — > U u (l) x Ud(l) X U s (l). 

• The EM interactions break all symmetry subgroups down to SUj(2) x U c (l) x Ub(l) x 
U U {1) x U d {\) x U S {1). 

The analysis in this work has the objective of sorting out these three sources of symmetry 
breaking from the current knowledge of the heavy meson masses. The analysis only involves 
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meson mass differences. Those in the D mesons are very well known with errors smaller than 

n 

2.5% [2j. In the B mesons, not all mass splittings are established; the mass differences (B*~ — 
B~) and (B*° — B°) are not separately known and, in addition, the errors are significantly 
larger than in the D-system |2|. Nonetheless, the available information is sufficient for the 
analysis to lead to significant conclusions. 

The analysis is based on the mass formuli that result from the expansions in 1 / itlq , where 
rriQ is the heavy quark mass (Q = c,b), in m u , and m s , and in the fine structure constant 
a. These mass formuli are similar to the ones given long ago by Rosner and Wise except 
that the QCD running of some of the parameters are taken into account. Using the notation 
H = M H , etc., they read as follows: 

ffM = Hi*\m Q ) + «W(fng) m q + a (a^(m Q ) Q 2 q - b^(m Q ) Q Q Q 9 ), (1) 

where q and Q denote respectively the light and heavy quark flavors of the heavy meson, 
the label * is used for the vector mesons, and Q g and Qq are the respectively the light and 
heavy quark charges. 

The first term in the mass formula is made out of the contributions in the limit m q — > 
and a — > plus SU(3) singlet contributions by light quark masses, i.e. given in terms of 
m u + ma + m s , and by electromagnetism. These SU(3) symmetric contributions are simply 
absorbed into the coefficients that determine //g*^(mg) in Eqn.(2) below. Our analysis 
cannot determine those contributions separately. The second term in the mass formula 
provides the SU (3) breaking effects by the quark masses to first order, and the last two terms 
provide the SU(3) breaking by the electromagnetic interaction, the first one representing the 
electromagnetic self-energy of the light antiquark, and the second one is the electromagnetic 
interaction between light antiquark and heavy quark. Throughout, mq represents the heavy 
quark pole mass. Up to (9(l/mg), the 1/mq expansion gives for Hq^ttiq): 

2 3 
TfL 771 

H*(m Q ) = H (m Q ) + h x (m Q ) + h 2 (m Q ) -f. (2) 

m Q m Q 

For dimensional purposes, we use the p-meson mass as the reference mass scale. The ex- 
plicitly displayed terms 0{1/w,q) and higher represent the mass difference produced by the 
strong hyperfine interaction. h\{rriQ) has a non-trivial dependence on mq that results from 
the QCD running of the effective heavy quark operator in the I/tuq expansion associated 
with it. The heavy quark effective Lagrangian in the 1/toq expansion is given in standard 
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notation by 



C eJf = Q v ivDQ v --l—Q v D 2 Q V + g s Q v G»» Q v + 0(l/m 2 Q ), (3) 

ZiriQ ^vtiq 

where the third term, which couples to the heavy quark spin, gives rise to h\. The coefficient 
Cm(p), where fj, is the renormalization scale, is given by j^]: 

C A 

c M (n) = c m Qm>) f^r41 00 (4) 



with the initial condition given by the matching to full QCD jjj: 

C M (m Q ) = 1 + (C A + C F ) ^4^- (5) 

Here Ca = 3,CV = 4/3, and (3q = 11 — |iV/ is the first coefficient of the /3-function with 
iVy the number of flavors lighter than irtq. Because h\ is proportional to Cm, we can express 
it in the following renormalization group invariant form: 

hi(m Q ) = C M {m Q ) a s (m Q ) «o hi, (6) 

where h% is rug and /i independent. In what follows, all over- lined coefficients are tuq and 
H independent. The term proportional to h 2 receives contributions from two terms in the 
heavy quark Lagrangian at 0(l/mn) [4]. Taking into account the running of h 2 with rriQ 
is thus impossible in this analysis. Fortunately, this is not important because the 0(l/mg) 
terms play a minor role in the B mesons, and thus neglecting the running of the h 2 is a good 
approximation. 



For the light quark mass effects, vn q {\i) is defined in MS scheme and one has: 

K(m Q ,n) = «o(^) + «i(/i) — - 7 /4(m Q ,/i) — , 

mg 4 m Q 

??7> 1 77Z 

K*(m Q ,fi) = «o(/x) + — + 7 Ki(m Q ,/x) — ^ . (7) 

m Q 4 m Q 

Because the spin-independent term in the (9(l/mg) heavy quark Lagrangian is scale inde- 
pendent, Ko and K\ are independent oirriQ and their dependence on is given by the running 
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of m 9 (/x), i.e. = ft> /q, i = 0, 1. On the other hand, because is proportional 

to Cm, it has an extra running factor similar to that of hi. 

K[(m Q , fj) = C M (m Q ) a 8 (m Q ) Po a a (fj) ^ R x ' . (8) 
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The fj, dependence given above is rather immaterial in our analysis, where the same // is 
used for D and B mesons. 

In the case of electromegnetic effects, we need to discuss seperately the two terms. The 
coefficients of the self-energy can be expressed in the following most general form: 

a(mn) = a + a 1 a^rriQ) 

TTlQ 4 ITiQ 

a*(m Q ) = a + ai — + 7 — ■ ( 9 ) 

m Q 4 ^ rriQ 

We note that a[ will run with uiq in a similar form as h±. The spin of the heavy quark will 
affect very little the self-energy term, and therefore a[ will be small. As shown later, a[ can 
be eliminated because of linear dependencies, which means that the effect associated with 
it cannot be distinguished from other effects on the meson masses. 

On the other hand, the electromagnetic terms involving the interaction between the light 
degrees of freedom and the heavy quark have the form: 

b(m Q )=b + b l ^- 3 -b[(m Q )^- 
ttiq 4 m Q 

b*(m Q ) = b + b, ^ + \ 6i(mg) ^ (10) 
m Q A m Q 

b and &i give the leading and subleading in l/mg contributions to the Coulomb interaction, 
and both are scale independent as it is known from the renormalization of current operators 
of the heavy quark, such as the electromagnetic current. The electromagnetic hyperfine 
effect proportional to b[ can be expressed as follows: 

b[(m Q , a a (n)) = b^ + C M {m Q ) a s (m Q )W b? . (11) 

It receives contributions from two general types of diagrams shown in Fig (1). The first 
term results from the coupling of the photon to the heavy quark spin and therefore it has no 
itlq dependence, while the second term corresponds to the coupling of a gluon to the heavy 
quark spin, and is therefore proportional to Cm and thus mg dependent. 

II. ANALYSIS 

In this analysis, we consider the five different mass splittings possible in each multiplet, 
i.e. (D + - D°), (D* + - D+), (D s - D°), (D*° - D°), (D* - D 8 ) and similarly for B-system. 
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The mass formuli leave one parameter independent mass relation, which reads: 

((B* s - B.) - (B*° - B )) = mcX . (m "| (p: - D.) - (D* + - D+)), (12) 

m b x{m c ) 

c A 

where we denote x{ m o) — Cm^q) OL s {mo) /3 ° • This mass relation is violated by terms 
0(1/i71q), 0(ot/mg), C(m^ 2 ), and 0(m g /mg). Note that this relation was discovered in 
Ref. [lj], where the evolution factor x{ m o) was n °t included. 

If one disregards the term 0(l/mn) in the strong hyperfme interaction, i.e. the term 
proportional to h,2, one obtains an additional relation: 

m cX (m b ) ((D*°-D°) + 2(D:-D S )) = 

m b X (m c ) {{B*- -B-)+2 (B* - B s )) ' 1 ) 

The deviations from this relation are a measure of the importance of the 1/rng term in the 
hyperfme interaction. 

In the mass formuli, there is a total of twelve parameters that enter in mass differences. 
Since there are ten mass differences and one parameter free mass relation, there must be three 
linearly dependent terms in the mass formuli that we can eliminate. The linear dependencies 
are such that di and b\ can be absorbed into do and bo, and a'^rriQ) into h\{rriQ) and b'^mq). 
Since no dependencies appear if one stays at leading order in 1/rriQ, it is natural to eliminate 
sub- leading terms. We, therefore, eliminate therefore di, a[ and b\. The linear dependencies 
imply that, in this analysis, one cannot determine the I/tuq corrections to the self-energy 
and to the Coulomb effects independently from other effects. In what follows, we therefore 
set: 

ai = ai(mg) = 6i = 0. (14) 

The quark masses are the key input parameters in the mass formuli. If we would only 
keep up to 0{1/w,q) terms, the ratio (D* — D)/(B* —B) would determine m c jm b = 0.40 for 
rrib ~ 5 GeV. There are more accurate determinations of the ratio m c / nib from the analysis 
of charmonium and bottomonium [5( that give m c /mb = 0.35 ± 0.03 and m& = 4.98 ± 0.13 
GeV. With this ratio, it becomes necessary to include the 0(l/mn) term in Eqn. (2). 
Note that Eqn. (12) could be used to extract the ratio m c /m b , but it requires a precision 
in the mass differences involved that is far beyond the current precision. For the light 
quark mass inputs, we only need their ratios, namely, m s /m where rh = (m u + md)/2, 
and m s /(md — m u ), which have been extensively studied in chiral perturbation theory j|. 
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The first ratio is obtained from the ratio M^/M^, which after next to leading order chiral 
corrections, gives a value of 24.4 ± 1.5, and the second ratio requires the input of isospin 
breaking observables, in particular the mass ratio (M K o — M K +)/M] c , and the rates for 
i] — > 3tt, with a result m s / (m^ — m u ) = 42.5 ± 3.2 [6|. For further reference, this latter ratio 
corresponds to having the electromagnetic mass difference (K + — K°)em — 2.0 ± 0.4 MeV. 
The final input is Aqcd required by a s ; we use Aqcd = 200 MeV. 

With our inputs for the heavy quark masses, we obtain for the left hand side of Eqn. (13) 
a value equal to 0.90 ±0.08, which gives some evidence for the need of the 0(1/171%) term in 
Eqn. (2). For Eqn. (13) to hold, it would be necessary to have m c /mb = 0.40. From Eqn. 
(12) we obtain the combination: 

(B*° - B°) - (B* - B s ) = -0.90 ± 0.16 MeV, (15) 

which is not known experimentally because the mass difference (B*° — B°) has not been 
established separately from the one for charged ones. The improvement over similar predic- 
tion given in l| is primarily due to the improved accuracy of the various inputs, especially 
the heavy quark masses and the running effect characterized by the factor x( m o)- 

In Table I, we give the results of our fits, displaying the partial contributions, and in 
Table II, the effects are combined. For the strong hyperfine contributions, we see that the 
0(1/i71q) effect tends to reduce the contribution from the leading O(l/mo) term by up to 
20% in D-mesons and by up to 8% in B-mesons. The errors quoted for the individual terms 
are rather large which indicate a strong correlation. This is displayed in Fig. (2). The error 
due to the uncertainty in m c /m{, manifests itself in the individual strong hyperfine terms 
where it is approximately ±30MeV in the D mesons. This indicates that, at present, the 
(9(l/mg) strong hyperfine effects cannot be established with good precision. Clearly, the 
combined hyperfine effects are very precise as shown in Table II. 

The leading order contributions by the light quark masses, which are independent of 
the meson's spin, i.e. the contributions proportional to k , are well determined by the 
fit, the error being dominated by the errors in the input light quark mass ratios. The 
0(1/ mo) corrections proportional to n\ can also be determined quite easily: because the spin 
dependent term k>2 turns out to be rather small, the K\ term is almost entirely determined 
in terms of the combination (D s — D + ) — (B s — B°). In particular, its sign is positive 
because (D s — D + ) > (B s — B°). This is a puzzling fact: if one tries to interpret the K\ 
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term from the point of view of a non-relativistic constituent quark model, the dependence 
of (H s — H) on the light quark masses will be through the reduced mass of the qQ system, 
and because the 0{m q /mQ) correction to the reduced mass is negative, the sign of the K\ 
term should be opposite to that of the k term. This is indication of a clear departure from 
the non-relativistic quark model picture which we have not seen addressed in the literature. 
To complete the light quark mass effects, the terms proportional to H2 provide the SU (3) 
breaking effects to the strong hyperfine interaction , and are determined with an error of 
about 10%, which is remarkable. For instance, in the case of the (D* — D s ) mass difference, 
it gives an upward shift of 3.3 ± 0.3 MeV. The combined effects of the light quark masses, 
displayed in Table II, show that they are determined by this analysis with an accuracy 
that is in general better than 10%. There is, however, one important and still unresolved 
problem concerning the light quark mass effects, and this has to do with the non-analytic 
contributions proportional to 0(m^ 2 ) 7], which are expected to be large sj], according to 
the estimated value of the coupling g [9J that gives the amplitudes D* — > Dir. A consistent 
analysis should include up to the m 2 q effects, which is beyond the current analysis. This 
problem, therefore, introduces some uncertainty in the determination of the light quark 
mass effects that is difficult to estimate. 

As mentioned earlier, we can separate the electromagnetic effects into self-energy and 
Coulomb plus hyperfine type terms. The self-energy has only a spin independent piece 
and it represents an effect of less than 1 MeV, and is determined with about 30% error. 
It has the same sign and comparable magnitude to results from calculations based on the 

nn 

Cottingham formula for electromagnetic mass shifts in a VMD approximation [10, HI[ • The 
effect of the Coulomb interaction is given by the term proportional to bo and its subleading 
piece proportional to bi, as explained earlier, has been absorbed into other terms. The fit 
determines the Coulomb effect with an error of 11%. The electromagnetic hyperfine effects 
are associated with the two parameters frj 1 and b^ . The input for (B* — B) does not 
differentiate between the neutral and charged mesons; if it is identified with the neutral 
mesons, then one of these two parameters can be eliminated as a consequence of the relation 
Eqn. (12). For this reason, in the fits, we carried out identifying (B* — B) = (B*~ — B~ + 
B*° — B°)/2; one finds a large correlation between these two parameters, which requires that 
we keep only one of them. We keep b[ , which amounts to ignoring the mg dependence 
in b^\ absorbing the rest of it into b[ . Our analysis is, therefore, insensitive to the QCD 
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running of the electromagnetic hyperfine effects, which is not surprising. We have moreover 
checked that our results are almost insensitive to the interpretation of the input (B* — B) 
as an arbitrary combination of the charged and neutral mass differences. The hyperfine 
effects are significant in the D-mesons, for instance, in the (D + — D°) case, it is about 60% 
of the Coulomb effect. On the other hand, in the B-system, the hyperfine effects are much 
smaller than the experimental uncertainties in the mass differences. Comparison, with the 



calculations in Ref. 



Ill ], shows agreement with the results in the elastic approximation to the 



Cottingham formula using VMD. Note that the inelastic contributions in the Cottingham 
formula that correspond to the interaction Coulomb and hyperfine terms are suppressed by 
1/tuq and, therefore, to the order we are working here, they can be neglected. 

It is instructive to make some comparisons. The electromagnetic shifts for the pseu- 
doscalar mesons are (D+ - D°) EM = 2.3 ± 0.2 MeV, and (B~ - B°) EM = 1.9 ± 0.2 MeV, 
which are similar within errors to (K~ — K°) E m — 2.0 ± 0.4 MeV. This is however a 
bit of a coincidence, as it can be seen from the results obtained for the different contri- 
butions, which are in the Coulomb and hyperfine cases very different in the D and B 
mesons. If we consider the vector mesons, we obtain (D* + — D*°) E m = 0.8 ± 0.2 MeV, 
and (B*- - B*°) EM = 1.6 ± 0.2 MeV, which we can compare with (p+ - p°) = 0.7 ± 0.8 
MeV (note that, as in the case of pions, the mass difference between charged and neutral 
p meson is purely electromagnetic). One important observation is that the uncertainties 
in the ratios of light quark masses only have a noticeable effect in the self-energy terms, 
while the uncertainty in m c /mf, leaves the results for the electromagnetic effects virtually 
unchanged. The isospin mass splittings due to m u and are: (D + — D°) mq = 2.42 ± 0.12 
MeV, (D* + -D*°) mq = 2.50±0.12 MeV, (B°-B~) mq = 2.16±0.2 MeV, and (B* -B*~) mq = 
2. 18 ±0.2 MeV. The difference between the last two is the negligible isospin breaking induced 
on the strong hyperfine interaction in the B mesons. This effect is larger and significant 
in the D mesons as shown by the difference between the first two mass differences. Fi- 
nally, one can give an accurate prediction: (B*° — B*~) = 0.6 ± 0.2 MeV. The combination, 
f (§(£>+ - D°) + (D* + - D*°) - l(B° - B~) - (B*° - B*-)) mq = 0.3 ± 0.02 MeV, gives the 
difference between D and B mesons of the O{l/mo) spin independent part of the isospin 
breaking by the quark masses. The electromagnetic part of the same combination is sub- 
stantially larger and equal to 2.8 ± 0.2 MeV. A similar analysis is straightforward for the 
mass splittings due to m s . 
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In summary, we have analyzed the different contributions to the mass splittings in heavy 
ground state mesons. The analysis shows that, with the current empirical accuracy of the 
heavy meson masses, one can determine these contributions with significant precision even 
at the sub- leading order in l/mq. The results obtained here can be useful for constraining 
models of heavy mesons, and perhaps also for lattice QCD calculations of heavy mesons 
masses where it is possible to study the light quark mass dependence. 
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AM 


hi 


h 2 




Kl 


< 


ao 


So 


B« 


D+ -D° 








2.01 ±0.14 


0.47 ±0.06 


-0.06 ± 0.01 


-0.7 ±0.2 


1.88 ±0.20 


1.17 ±0.06 


D s -D+ 








82.2 ±2.8 


19 ± 2.3 


-2.3 ±0.2 











D*° - D° 


175 ± 29 


-34 ± 28 








0.09 ±0.01 








1.0 ±0.05 


D*+ -D+ 


175 ± 29 


-34 ± 28 








0.17 ±0.01 








-0.52 ±0.03 


D* s - D a 


175 ± 29 


-34 ± 28 








3.3 ±0.28 








-0.52 ±0.03 


B° - B~ 








2.01 ±0.14 


0.16 ±0.03 


-0.016 ± 0.002 


-0.7 ±0.25 


-0.94 ±0.1 


-0.20 ±0.02 


B* - B 


49.8 ±3.5 


-4.1 ± 3.3 








0.038 ± 0.005 








-0.05 ±0.01 


B s -B 








83.2 ±2.9 


6.8 ± 1.2 


-0.68 ± 0.08 


-0.36 ±0.1 


-0.47 ±0.05 


-0.10 ±0.01 


B* - B a 


49.8 ±3.3 


-4.1 ± 3.3 








0.94 ±0.13 








0.09 ± 0.01 



TABLE I: The mass shifts due to the different terms are labeled by the corresponding coefficient 
in the mass formuli and are in the in units of MeV. The errors include the uncertainties in the 
input quark mass ratios. In the fit (B* — B) is interpreted as \{B*~ — B~ + B*° — B°). 



AM 


Strong HF 


Light quark masses 


Electromagnetic 


Total 


PDG [2] 


D+ - D° 





2.42 ±0.12 


2.33 ± 0.22 


4.75 ±0.15 


4.78 ±0.10 


D S ~D+ 





98.96 ± 0.49 





98.96 ±0.47 


98.85 ±0.30 


D*° - D° 


140.99 ±0.1 


0.09 ±0.01 


1.04 ±0.05 


142.12 ±0.1 


142.12 ±0.07 


D*+ - D+ 


140.99 ±0.1 


0.17±0.01 


-0.52 ±0.03 


140.64 ±0.09 


140.64 ±0.10 


D* - D a 


140.99 ±0.1 


3.30 ±0.28 


-0.52 ±0.03 


143.77 ±0.3 


143.8 ±0.4 


B° - B~ 





2.16 ±0.12 


-1.86 ±0.19 


0.304 ±0.11 


0.33 ±0.28 


B* - B 


45.70 ±0.02 


0.04 ± 0.004 


-0.05 ±0.002 


45.69 ± 0.02 


45.78 ±0.35 


Bs-B 





89.32 ±0.37 


-0.93 ±0.09 


88.39 ±0.33 


88.3 ±1.8 


B* - B s 


45.70 ± 0.02 


0.94 ±0.08 


0.09 ± 0.005 


46.73 ± 0.08 


45.3 ± 1.5 



TABLE II: Mass contributions by strong hyperfine, light quark masses and electromagnetism in 
units of MeV. The errors include the uncertainties in the quark mass ratios. The fit has % 2 ~ 1. 
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(a) 



(b) 



FIG. 1: Diagrams contributing to the electromagnetic mass shift parameters: type (a) contributes 
to , and type (b) to . The crosses indicate respectively the coupling of the photon and the 
gluon to the heavy quark spin. 




FIG. 2: Statistical correlation between the leading and sub-leading hyperfine contributions in D 
mesons for the fit with m c /mb = 0.35. 
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